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We theoretically study subgap states appearing at the interface between two three-dimensional 
topological insulators which have different configurations in the spin-orbit interactions from each 
other. The coupling of spin a with momenta p is configured by a material dependent 3x3 matrix 
A as a^A^piy. We show that the spectra of the interface suggap states depend strongly on the 
relative choices of A in the two topological insulators. In particular, we focus on properties of 
gapless states which appear when A in two topological insulators are connected by the inversion in 
momentum space. We also discuss the robustness of the gapless states under perturbations breaking 
the time-reversal symmetry or the band-inversion symmetry by the numerical simulation. 

PACS numbers: 73.20.At, 74.45.+C 



I. INTRODUCTION 

Topological insulators (TIs) are new class of condensed 
matter— .i. A topological number Z2 defined in terms of 
the global property of wave function for the occupied 
states below the gap distinguishes the topological insu- 
lating phase Z2 = 1 from the conventional one Z2 = 0. 
The bulk-edge correspondence guarantees the presence of 
gapless states at the surface of Tl^^^^. On the surface 
of three-dimensional TIs, the excitation spectra of the 
gapless state are described by so called Dirac cone, (i.e., 
E = The strong spin-orbit interaction locks the 

direction of momenta and that of spin, and is responsible 
for the non- trivial Z2 topological number. 

When we focus on the surface state of a single TI, the 
spectra of the surface state are independent of the con- 
figurations in the momentum-spin locking. On the other 
hand, when we focus on the interface states between two 
different topological insulators, the spectra of the inter- 
face state depend on the relative configuration of the 
spin-orbit coupling in the two TIs^^^. The investigation 
in this direction has been also done in superconducting 
junctions consisting of two helical superconductor^ and 
of two chiral superconductors^iii. These studies focus 
on the discrete degree of freedom such as helicity, chi- 
rality and mirror symmetryi^i^. However, the coupling 
of spin a with momenta p is configured by a material 
dependent 3x3 matrix A as a^A^pj^. The properties of 
the interface subgap states depend on the relative choices 
of A in the two TIs. Generally speaking, A in one TI is 
transformable to A in the other not only by the inver- 
sion but also by the rotation. The former is a discrete 
transformation but the latter is a continuous one. 

In this paper, we discuss the properties of the two- 
dimensional states appearing at the interface of two 
three-dimensional TIs which have different configura- 
tion in the spin-orbit coupling. The configuration of 
the spin-orbit coupling is mainly characterized by the 
two transformations: (i) the sign inversion in the spin- 
orbit coupling term for the momentum perpendicular to 
the interface plane and (ii) the continuous rotation of 



the momentum-spin locking angle {(j)) within the par- 
allel plane to the interface. In the analytical calcula- 
tion, we obtain the spectra of the interface subgap states 
in the presence of both the time-reversal symmetry and 
band-inversion symmetry. When the two topological in- 
sulators have the same sign in the transformation (i), 
the zero-energy surface states appear only at the unique 
angle of in the transformation (ii). The zero-energy 
states in this case are characterized by the mirror Chern 
number- and fragile under the continuous transforma- 
tion (ii). When the two topological insulators have the 
opposite sign to each other in the transformation (i), on 
the other hand, the gapless surface state appears for all 
angles in the transformation (ii). The gapless states for 
the latter case are protected by the Sato's winding num- 
ber—— and robust under the the continuous transfor- 
mation (ii). We also carry out numerical simulation on 
the three-dimensional tight-binding model to check the 
robustness of interface states at the zero-energy under 
the perturbations breaking the time-reversal symmetry 
or band-inversion symmetry. It has been already known 
that the surface state of a single TI are sensitive to the 
direction of Zeeman fieldi^ — . We find that the inter- 
face states at the zero-energy are also sensitive to the 
direction of Zeeman field. In the presence of the Zeeman 
field parallel to the interface plane, the zero-energy inter- 
face states protected by the mirror Chern number van- 
ish, whereas those protected by the Sato's winding num- 
ber remain unchanged. In the presence of the Zeeman 
field perpendicular to the interface plane, on the other 
hand, the zero-energy interface states protected by the 
mirror Chern number remain unchanged, whereas those 
protected by the Sato's winding number vanish. The per- 
turbations breaking the band inversion symmetry remove 
the interface zero-energy state in both cases. 

This paper is organized as follows. In Sec. II, we 
explain the theoretical model considered in this paper. 
The analytical expression for the spectra of the interface 
states is given. In Sec. Ill, we discuss analytical results 
of the interface states for two types of junctions. At the 
same time, we check the robustness of the zero-energy 
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interface states by the numerical simulation on the tight- 
binding model. The conclusion is given in Sec. IV. 



II. MODEL 

The most simple Hamiltonian of three-dimensional 
topological insulator reads, 

_({m-hp^)(j^ aap \ 
^^^-[ aap -{m-bp')ay^ 

where a, 6, and m are positive constants, M = {m — hp^)^ 
cr^ is 2 X 2 identity matrix, and cr = (a^^a^^a^) are 
Pauli matrices in spin space. The e igen values o f the 
Hamiltonian are itE'p with Ep = ^/M^^^^aP\p\^ . For 
convenience, we utilize the short notation as, 



(2) 



where M is regarded as the Dirac mass and we use 4x4 
Dirac matrices. 











-a^ 



Eq. ([2]) is called Dirac Hamiltonian and describes topo- 
logical materials such as topological superconductors and 
superfluid ^He-B phased—. In real topological insu- 
lators, the coupling between spin and orbital parts has 
more general form 



Ho = aa^K^p^, + M/3, 



(3) 



where the configuration matrix defines the angle 
which locks spins and momenta and describes the rota- 
tion or the inversion in momentum space. In such case, 
K^^ is the real symmetric matrix satisfying AJ;^A^ = 3^. 
The Z2 topological number does not distinct two topo- 
logical insulators with different A^^ because they are con- 
nected with each other by a unitary transformation^. 
Thus the energy spectra of Eq. (|3j) are independent of 
A^. 

When we focus only on a isolated topological insulator, 
physics of of Eq. (j2j) is the same as that of the simple 
Hamiltonian in Eq. (|2]). However, when we consider a 
junction of two different topological insulators, physics 
happening near the junction interface depends on the 
choice of A^^ in the two topological insulators. This is 
because there is no unitary transformation which trans- 
forms the two different K^^ into the Hamiltonian in Eq. [3] 
at the same time. In the following, we study the prop- 
erties of subgap states at the junction interface of two 
different topological insulators. We apply two theoretical 
methods to analyse the spectra of the interface states. At 
first, we analytically derive the energy spectra of the in- 
terface states from the boundary condition of wave func- 
tions. This method gives analytical expression of the 
subgap spectra of the interface states. Its applicability. 



however, is limited to the Hamiltonian preserving of the 
band-inversion and the time-reversal symmetry. In addi- 
tion to the analytical method, we also perform numerical 
simulation on the three-dimensional tight-binding model 
to check the robustness of the predicted gapless states by 
the analytic method. They are complementary methods 
with each other and gives the consistent results. 



A. Spectra of interfacial bound states 

We begin with the Hamiltonian in Eq. (|3]) which is 
represented as 

Mcr^ aa^^A^px 



Ho 



(4) 



The Hamiltonian in Eq. (j4]) preserves the band-inversion 
symmetry. 



VH{p)V-^ =-H''{-p), 
V 



-(jy 
(jy 



The wave function can be given by 



*(+)(r) 



-VpV 




(5) 
(6) 

\ (7) 
(8) 



where ^(+) and are the wave functions belonging 

to Ep and —Ep^ respectively. It is evident that Ep is 
independent of A^. 

When two topological insulators touch with each other 
at z = 0, the wave function for ^ = E'p > are repre- 
sented by 



X e 



for 2; < and 

X e 

for 2; > with 



va 
uT 



(1) 



ihp 



ihp 



(10) 



(11) 



^ti) = T-yiH))tip±)>^^ p± = (fc,±fc.), (12) 



(13) 
(14) 
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p={k,k^), k = {kx,ky), p={x,y), (15) 

bp, (16) 



a" 
462' 



(17) 



where / = 1 and 2 indicate the two topological insulators. 
The 2x1 coefficient matrices A, 5, C, and D represent 
the amplitude of outgoing waves from the interface. At 
a surface of a topological insulator, it is easy to confirm 
that the spectra of the surface state become E = a\k\ 
from the boundary condition ^^^\r)\ 



0. 



At the interface of the two topological insulators, the 
boundary conditions at z = 0, 



gives a condition 



(18) 
(19) 



det 



-)W (r+)-r-)) (r^,)-r+)) 



(1) 



(2) 



0. 

(20) 



In the quasiclassical approximation, we use a relation 
<^ mb^ which allows — <C k^k~ ^ k^. Thus 



3 find 
det 



(1) 



(2) 



det 



77 



^ ^1) 



0. (21) 



From the condition, we obtain the energy of interface 
subgap states 



E = apF 



\ 



2|p|2 



(22) 



The results indicate the spectra of subgap interface states 
depend on relative choice of A^i) and A(2)- Within 
Eq. (|4]), the spectra obtained from Eq. ([2T]) shown in 
Eq. (|22]) always doubly degenerate. 



B. Numerical simulation on tight-binding model 

To confirm obtained results in the analytical method 
and to check the robustness of the interface gapless states 
under perturbations, we also perform the numerical cal- 
culation on the tight-binding model. We describe the 
topological insulator by using the two-band model as 



(23) 



M = (m + 2^2 (cos (/ca^ Co) + cos{kyCo)) — 2))6jjf 

-2M,y +^i((^,,y+i+(^,,/-i), (24) 
hso = 0.20-'' [A^ sin(/c^co) + A^ sm{kyCo)] Sjj^ 
-mia^A^((5,-,v+i -(5,-^.,_i), (25) 



where k = {kx^ky), and j represents the place in the 
2;-direction. We apply the periodic boundary condition 
in the X7/-plane and the hard wall boundary condition in 
the 2;-axis. We choose ai = 7.86m/co, a2 = 14.6m/co, 
bi = 3.57 X lOm/cg, 62 = 2.02 x lO^m/cg with the lattice 
constant cq being 5 [A] in both sides taking into account 
band structures of Bi2Se3^^. In the simulation, the length 
of the each topological insulator in the z-direction is 400 

Co. 

To check the robustness of the interfacial zero-energy 
states, we also consider two types of perturbations in 
addition to Eq. (|23l) 



^Bi(j,/) = 1 



4x4 



M(2)^j,j' + ^(2)(^j,j'+i + for j > 



^TR 



h • a 







h • a 



(26) 
(27) 



Eq. ([26|) represents the constant shift of the chemical po- 
tential and the band asymmetry in the two TI which 
break the band-inversion symmetry. Eq. (|27|) represents 
the uniform Zeeman potential which breaks the time- 
reversal symmetry. 



III. RESULTS 

In what follows, we fix A (2) at diag( 1,1,1) and study the 
spectra of the interface for several choices of A(i) because 
the relative configuration of momentum-spin locking is 
responsible for this physics. 



A. Rotation within the x^y-plane 



When we choose 
A(i) = 



COS ( 

sin^ 




— sin (f) 
cos(j) 
1 



the spectra in Eq. (|22]) become 



E = a 



\k\^sm^{^/2). 



(28) 



(29) 



The spin-orbit interaction in Eq. ([28]) describes only 
the continuous rotation of spin-momenta locking angle 
(f) within the xy plane. The spin- moment a locking in the 
z axis is common in the two topological insulators. In the 
effective model around the T point in the Brillouin zone, 
the zero-energy state is possible only at \k\ = pr and 
= TT. In Fig. [H we show numerical results of the dis- 
persion relations for several choices of continuous trans- 
formation of (/). The dispersions are plotted as a function 
of kx because the numerical results are totally isotropic 
in momentum space. It is evident that the interface sub- 
gap state is absent at = because the resulting relation 
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FIG. 1. The energy dispersion of interface states along kx 
is shown for Eq. ([28]). The dispersion has the rotational 
symmetry in two-dimensional momentum. We choose (f) 
0. 7r/2,and tt in (a), (b), and (c), respectively. In agreement 
with the expression in Eq. ([29]) . the zero-energy states ap- 
pears only = TT as shown in (c). 



distinguishes the two topological insulators and protects 
the zero-energy stated. However, our numerical calcu- 
lation suggests that the gapless state appears only at 
(j) = TT. Therefore the delicate material tuning is nec- 
essary for having the ring-shaped zero-energy state. A 
small perturbation modifying the locking angle (f) and 
may open the gap. In fact, Takahashi et al.^ has tried to 
find reasonable combinations of two different topological 
insulators such as Bi2Se3, Bi2Te3, and Sb2Te3. 

Next we check the robustness of the gapless inter- 
face state at ^ = TT by introducing the Zeeman effect 
in Eq. (|27|) which breaks time reversal symmetry. In 




FIG. 2. The energy dispersion of interface states along kx 
under the Zeeman field are shown for Eq. ()28p with cj) — 
The figure (a) and (b) are calculated in the presence of the 
Zeeman field of hx — 0.5m and hz — 0.5m, respectively. 



A(2) = A(i) means the junction of two identical topolog- 
ical insulators. The gap of the interface state decreases 
with increasing (/) to tt as shown in Figs. 1(a) and |l(b) 



Fig. 2(a)[ we consider the Zeeman field parallel to the 



The wave numbers at the minima in the upper band and 
the maxima in the lower band have the ring-shaped in two 
dimensional Brillouin zone. At ^ = tt, such minima and 
maxima are touched with each other as shown in Fig l(c)[ 
The surface state becomes gapless and has ring-shaped 
Fermi surface. The dispersion along radial momentum 
represent a linear dispersion from the ring-shaped zeros. 
Such ring-shaped zeros have been discussed in Ref. 0. 
The choice of ^ = tt corresponds to the inversion of the 
two axes within the interface plane, which results in the 
opposite helicities in the two topological insulators. A 
topological number derived from the Mirror symmetry 



interface plane with hx = 0.5m and hy = hz =0. The 
interface state between the two topological insulators be- 
comes gapped in the presence of the in-plane Zeeman 
field. The Zeeman field perpendicular to the x^-plane, 
on the other hand, dose not remove the gapless inter- 
face state shown in Fig. |2(b)[ where hz = 0.5m and 
hx = hy = 0. Thus the robustness of the interface zero- 
energy state depends on the direction of the Zeeman field. 
We note that the tendency in the magnetic anisotropic 
effect of the interface states is opposite to that found in 
the gapless surface state of a single topological insula- 



tor' 



21.27 



This discrepancy, however, is not surprising at 
all. The helicity is defined in terms of the spin-orbit cou- 
pling including the in-plane momenta in Eq. (j3j). Thus 
the perpendicular Zeeman field does not affect the spin- 
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orbit coupling between parallel momenta and spin. 




FIG. 3. The energy dispersion of interface states along kx 
with magnetization is shown for Eq. (|28p with cj) — The 
figure are calculated in the presence of the band asymmetry 
as — ~M(2) — 0.1m and — O.l^i. 



O.I61 



In Fig. 3(a), we calculate the energy dispersion at 
TT in the presence of the band asymmetry in Eq. 
where we introduce = — /i(2) = 0.1m, 6(1) 
and 6(2) = in Eq. (|3]). Obviously, the band asymmetry 
removes the gapless interface state. Therefore, it would 
be difficult to find the gapless state at = tt in Eq. ([28|) 
in realistic junctions. 



B. Inversion of the 2;-axis 



Next we choose 



(1) 



cos ( 
sin ^ 




— sm ( 
cos (j) 







-1 



(30) 



The spin-orbit interaction in Eq. ([30]) describes the in- 
version in the z axis plus the continuous rotation of spin- 
momenta locking angle (j) within the xy plane. The inver- 
sion of the z axis cannot be described by the continuous 
transformation. The spectra in Eq. (|22]) becomes 



(31) 



E = a|A^||cos((/)/2)|. 



At (/) 7^ TT, the interface state is described by the Dirac 
cone with the Dirac node at = 0. We note that this 
conclusion holds for all (f). Thus realizing the Dirac cone 
like interface gapless state in Eq. (|3Q|) is rather easier than 
that of the ring-shaped zeros in Eq. (|28|) with (j) = tt. 

At /c = 0, Eq. ([4]) is block diagonal for each spin space. 
It is possible to define the Sato's one-dimensional wind- 
ing number in each spin sector—"—. Details are shown 
in Appendix. By using the winding number defined in 
each spin sector, it is possible to topologically distinct 
the two topological insulators. The topological numbers 
are calculated as 




FIG. 4. The energy dispersion of interface states along kx is 
shown for Eq. (|3Q|) . We choose = 0, 7r/2,and tt in (a), (b), 
and (c), respectively. The gapless interface state appears for 
all 0, which shows the robustness of the gapless state against 
the changing the rotation angle (f). 



for z < and 



(2),t 



1, w, 



(2),i 



-1 



(33) 



for z > 0. The bulk-boundary correspondence suggests 
the number of zero-energy states for spin-up space is 



W(i),t = -1, W(i),i = l 



equal to — W(2),tl — ^ ^^^^ ^^r spin-down one 

is |W(i),| — W(2),il =2. As a consequence, four states 
are degenerate at the zero-energy and at /c = 0. 

The analytical results are confirmed by numerical cal- 
culation in Fig. m where we plot the dispersions as a 
function of kx for Eq. (|3Q|) with (j) = 0,7r/2 and tt. The 
numerical results in Fig.HJ^a) show the doubly degenerate 
Dirac cones at the interface in argument with the ana- 
lytical results. In Fig.HJ^b), the numerical results deviate 
from the analytical one and the two Dirac cones separate 
from each other as \kx\ increases from zero. The ana- 
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lytical results in Eq. (jSTJ) indicate the fiat zero-energy 
interface states at = tt for all k. However, the quasi- 
classical approximation is not valid at = tt. In fact, the 
two determinants of matrices in Eq. ([2T]) are automati- 
cally proportional to at = tt. In such case, we need 
to analyse more general condition in Eq. ([20]) beyond 
the quasiclassical approximation. The numerical results 
in Fig. 2(c) also show the existence of four zero-energy 
state at ^ = TT. One of the branch shows the linear dis- 
persion, whereas the other seems to have the quadratic 
dispersion from A; = as shown in Fig. 4(c). We note 
that the choice of = tt corresponds to the inversion in 
all axes. 




FIG. 5. The energy dispersion of interface states along kx 
with magnetization is shown for Eq. ()30p at = 0. The 
figure (a) and (b) are calculated in the presence of Zeeman 
field in the xy plane hx = 0.5m and that in the perpendicular 
to the xy plane hz = 0.5m, respectively. 

We check the robustness of the zero-energy states un- 
der the perturbations breaking the time reversal symme- 
try in Eq. (|27|) . In Fig. 5(a) and (b), we show the disper- 
sion of the interface states in the presence of the Zeeman 
fields within the xy plane with = O.bm^hy = and 
hz = in (a) and the Zeeman field perpendicular to the 
xy plane with hx = hy = and hz = 0.5m in (b). The 
results in Fig. 5(a)| show that the two Dirac nodes stay 
at the zero-energy and at /c = in the Brillouin zone 
even in the in-plane Zeeman field. The Sato's number 
cannot be well defined any longer in the absence of the 
time-reversal symmetry. However, the zero-energy states 
under Eq. ([30]) are very robust even in the in-plane Zee- 
man field. When we apply the Zeeman field along the z- 
axis, on the other hand, the interface zero-energy states 



5(b)[ Therefore the robustness of 



vanish as shown in Fig. 

the zero-energy states also depends on the direction of the 
Zeeman field. The zero-energy states are sensitive (insen- 
sitive) to the Zeeman field in the perpendicular (parallel) 
Zeeman field to the xy plane, which is the opposite ten- 
dency to that in Fig. 2(a) and (b). But such tendency 
is consistent with that found in the topological gapless 
states at the surface of ^He-B phase^^^^. Indeed, by the 
unitary transformation, Eq. (j4]) is deformed into 



/7th 



1 

-iay 

-ia^aa • p 



(34) 



(35) 



The last Hamiltonian is nothing other than the Hamil- 
tonian of ^He-B phase^^. The surface Majorana state 
of ^He-B phase is protected not only by Z2 number but 
also by the Sato's one. The common magnetic anisotropy 
found in the two systems might be originated from the 
same topological nature of the wave function. 




FIG. 6. The energy dispersion of interface states along kx is 
shown for Eq. (|3Q|) at = 0. The results are calculated in the 
presence of chemical potentials of = ~M(2) = 0.2m and 



the band asymmetry bm — O.lbi and 



0(2) 



0. 



In Fig. [6l we discuss effects of the perturbation break- 
ing the band-inversion symmetry in Eq. (|26|) with /i^i) = 
~M(2) = —0.2m, 6(1) = O.I61 and 6(2) = 0- The two Dirac 
cones are lifted off from each other and the gap is open- 
ing in the absence of the band-inversion symmetry. The 
Sato's number cannot be defined in the absence of the 
band inversion symmetry. Therefore there is no reason 
for the interface states remaining gapless. 



C. Remaining configurations 

Finally we choose 



A 



(1) 



s 
cos( 
sine 





— sin ( 
cos 6 



(36) 



with s = ±1. For s = 1, the transformation represents 
only the continuous rotation within yz-pldiiie. The gap- 
less states appear only at (j) = ir. Because x and y axes 
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are no longer equivalent to each other, the dispersion be- 
comes anisotropic in the two-dimensional Brillouin zone 
as shown in Figs. 7(a)| and |7(b)[ Here we plo t the disper- 
sion as a function of kx at ky = in Fig. |7(a) The results 
suggest the doubly degenerate Dirac cones. On the other 



hand, Fig. 7(a) shows the dispersion along ky at kx = 0. 



The two dispersion branches clearly separate from each 
other. The Dirac nodes, however, stay at A; = 0. 




kx[co ^] 




FIG. 7. The energy dispersion of interface states in Eq. ()36|) 
with s = 1 and (j) — In (a), /cy is fixed at the origin and 
the energy is plotted as a function of kx. In (b), kx is fixed 
at the origin and the energy is plotted as a function of ky. 



When we choose 5 = — 1, Eq. (|36|) represents the in- 
version in the x axis. The choice of = tt in Eq. (|36]) is 
identical to Eq. ([3Q|) with (/) = tt whose results are shown 
in Fig. |4(c) Therefore we seek gapless states for (/> ^ tt. 
The results for = are shown in Figs. |8(a)| and |8(b)[ 
We p lot th e dispersion as a function of kx at /c^y = in 
Figs. 8(a)[ The results suggest that there are the two 
Dirac cones but their nodes stay at finite value of kx- 
Correspondingly, there is gapless state in the dispersion 
along t he ky at the Dirac point of kx = ko = 0.0947r 
The similar dispersion can be seen also at 
At s 



in Fig. [8(b)_ 

(j) = i/T/Tls shown in Fig. [8(c)[ At s = -1, Eq. (j36 
includes the inversion in x direction. In such case, two 
Dirac cones appear with their nodes staying on kx axis 
and these cones come close to the origin of /c = as 
increasing of to tt in Fig. 4(c)[ 

The rotation in xz-plane is equivalent to that in yz- 
plane under interchanging kx and ky. Thus the Dirac 
nodes are stay on the ky axis when we consider inversion 
in the y direction. Together with the results in Fig. 4(c)[ 



we conclude that the gapless states appear when the rel- 
ative configuration matrix A(i) includes the inversion in 
odd number axes. Unfortunately, these gapless states are 
also fragile under the perturbation which breaks the-band 
inversion symmetry. 




7r/2 





FIG. 8. The energy dispersion of interface states along (a) kx 
and (b) ky are shown at = and s = — 1. In the figure (b), 
we set kx at the Dirac point ko — 0.0947r. The figure (c) is 
energy dispersion along kx at (j) = n /2 and s = 1. 



IV. CONCLUSION 

We have studied the interface state between two topo- 
logical insulators with different configurations of spin- 
orbit interactions. The two topological insulators are 
touching at 2: = and its interface is flat within the xy 
plane. The coupling of spin cr with momenta p is config- 
ured by a material dependent 3x3 matrix A as a^A^pi^. 
The spectra of subgap states depends on relative config- 
uration matrix A(i) = A(i)A^^, where A(/) for / = 1 and 
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2 characterize the spin-orbit couphng in the two topolog- 
ical insulators. The properties of interface gapless states 
in the presence of the time-reversal symmetry and the 
band- inversion symmetry are summarized as follows. 

(i) When A12 includes the inversion along one-axis, 
two Dirac cones appears in the two-dimensional Brillouin 
zone. When the inversion axis is in z direction, the dou- 
bly degenerate Dirac cones appear and their nodes stay 
at the F-point (kr^. = ky = 0). When the inversion axis is 
in X {y) direction, the two Dirac cones are separated in 
the Brillouin zone and their nodes stay on x {y) axis. In 
this case, the Dirac nodes are robust even if A12 includes 
the rotation around the inversion axis. 

(ii) When A12 includes the inversion in two-axes, the 
spectra depends also on the choices of the two-axes. 
When the two- axes are x and the junction have the 
mirror symmetry and the ring-shaped zeros appear in 
the Brillouin zone. When one of the two- axes is z and 
the other is x (?/), two Dirac nodes appear at the F-point. 
The dispersion of them are degenerate along ky {kx) but 
the degeneracy is hfted along kx {ky). In this case, the 
Dirac nodes are fragile under the rotation A12 around 
any axis. 

(iii) When A12 includes the inversion in all axes, non- 
degenerate two Dirac cones appear and their nodes stay 
at the F-point. This situation can be regarded as the 
series transformation of (i) and (ii). Therefore, the robust 
gapless state against rotation of (j) is realised under the 
inversion in odd number of axes. 

The numerical results show that the interface gapless 
states are sensitive to the direction of the Zeeman field 
breaking the time-reversal symmetry. They are robust 
under the Zeeman field in one direction, but fragile un- 
der the Zeeman field in another direction. Unfortunately, 
all gapless states are fragile in the presence of the per- 
turbations breaking the band-inversion symmetry. The 
last property implies the difficulty of finding the gapless 
states within the combination of existing topological in- 
sulators. However, our results predict unusual gapless 
state appearing at the interface of two different topolog- 
ical superconductors and that of two different superfluid 
phases. This is because Bogoliubov-de Gennes Hamilto- 
nian similar to Eq. ([3j) always satisfies the particle-hole 
symmetry (band-inversion symm etry in this paper) so 
called as chiral symmetry in Ref. 
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Appendix A: Sato's winding number 



When we choose 



A 



COS( 

sine; 




— sm ( 

COS(j) 





(Al) 



with s = ±1, the Hamiltonian in Eq. (j4]) at /c = is block 
diagonal in each spin space irrespective of (j). Namely 



askz 



askz 



Hi 



-askz 



-askz 
-Mz 



, (A2) 



with Mz = m — bk^. For instance, eigen values of are 
±e/e^ and eigen vectors are represented by 



and 



for and — e/c^, respectively. Here we define 



cos Uk 



(A3) 

(A4) 
(A5) 



In the presence of the time-reversal symmetry, the eigen 
vectors can be represented only by real quantities. In 
the presence of the band-inversion symmetry, the Sato's 
winding number can be defined at /c = as 

W(t,5) = ^ f dkzdkjk.. (A6) 

To estimate the topological number, we use tight-binding 
representation of the Hamiltonian, 

Mz =m-hkl^m- 2t{l - cos kz), (A7) 
akz aot sin kz^ (A8) 

with t>0, m — 4t<0, and |ao| <C 1 is a dimensionless 
constant. We find that 



>V(t, s) = sgn(s), W(;, s) = -sgn(s). 



(A9) 
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